Design, Synthesis, Molecular Modeling, In Silico ADME Studies and Anti-HIV-1 Assay of New Diazocoumarin Derivatives.

Some new diazo incorporated coumarin compounds were designed and synthesized to evaluate their anti-HIV activity. Overall, compounds were active against HIV at 100 μM. Additionally, no cytotoxic effect was observed at this concentration. The compound with 4-chlorobenzyl group indicated the best anti-HIV activity (52%). Docking studies using the later crystallographic data available for PFV integrase showed similar binding modes to HIV-1 integrase inhibitors. On the basis of these data, nitrogen atoms of 1,3,4-oxadiazole ring have been involved in the Mg2+ chelation and 4-chlorobenzyl group occupies the same position as 4-flourobenzyl group of raltegravir in the active site. In addition, in silico ADME assay demonstrated favorable physicochemical properties for the new designed compounds. Thus, synthesized structures could be introduced as a novel template for designing safe anti-HIV compounds with integrase inhibitory potential.


Introduction
Acquired Immunodeficiency Syndrome (AIDS) is the second cause of death from infectious diseases in low income countries (1). Human Immunodeficiency Virus (HIV) is recognized as the infective agent. Current anti-HIV drugs target different stages of HIV life cycle including absorption, fusion, reverse transcription, integration, and proteolytic cleavage (2-4). Although there is a plethora of anti-HIV drugs, drug toxicity, viral resistance, and patient incompliance has been limited anti-HIV therapy (5). Altogether, there is a demand for novel drugs from alternative scaffolds. HIV-1 integrase is an essential enzyme which catalyzes the integration of viral DNA into the host cellular DNA. This process is separated into 2 steps. First, 3′-processing resulting in the creation of free hydroxyl groups at reversetranscribed proviral cDNA 3′-ends. During the second step (strand transfer), integrase catalyzes the insertion of the 3′-processed viral DNA into the host cell genome (6, 7). FDAapproved integrase inhibitors (raltegravir (8, 9), elvitegravir (10, 11) and dolutegravir (12,13) in Figure 1) compete with host DNA in the IN active site, thus specifically inhibiting the strand transfer step (14). Development of integrase strand transfer inhibitors (INSTIs) represents an attractive research area in the field of HIV-1 drug discovery because there is no host cell equivalent for HIV-1 integrase enzyme (15).
The integrase strand transfer inhibitors (INSTIs) bind in the interface of the viral DNA and catalytic DDE motif. For this reason, they could be classified as interfacial inhibitors. The DDE motif contains a triad of three conserved acidic residues (Asp-64,  coordinating 2 Mg 2+ cofactors which are critical for the catalytic activity (16). There are various chemical classes of INSTIs such as oligonucleotides, peptides, hydroxylated aromatics, hydrazides, quinoline, and quinolones. All specific INSTIs feature a metal chelating motif that complexes 2 Mg 2+ cofactors and a hydrophobic terminal benzyl moiety (17)(18)(19).
Coumarin is a well-known drug scaffold which displays versatile biological activities especially substantial anti-HIV activity (20). Notably, IN inhibitory activity has been reported for coumarin dimmers (1, Figure 2) (21) and 4-hydroxy-5-azacoumarin-3-carbox(thio)amides (2, Figure 2) (22). Furthermore, the anti-HIV activity of a bis-azo compound (FP-21399, Figure 2) has been evidenced (23). With regard to these data we selected coumarin scaffold as a core moiety to design new anti-HIV compounds. The decision for selecting substitutions on coumarin ring was carried out from the information obtained through our previous studies. Recently, we reported development of 8-phenylquinolin-4(1H)-one derivatives possessing 1,3,4-oxadiazole ring as metal-chelating moiety along with a benzyl part necessary for HIV-1 IN inhibitory activity. Among them, compound 3 ( Figure 2) displayed promising anti-HIV-1 activity with EC 50 = 50 µM and no considerable cytotoxicity (CC 50 > 100 µM) (24). SAR analysis showed that the incorporation of electronegative substitutions like halogens (-F, -Cl) at proper positions of benzyl ring would improve inhibitory activity (25). In this study, we merged coumarin core with a 5-(halo-substituted) benzylthio-1,3,4oxadiazole motif through a diazo-phenyl linker (Designed structure, Figure 2) to present new structure for anti-HIV activity. Some compounds containing alkyl substitution such as methyl, ethyl, and propyl instead of benzyl ring at C-5 position of 1,3,4-oxadiazol ring were also designed. Using wide range of substitutions would help us to perform a SAR analysis on this series.
The overall shape and functionalities of designed structure particularly resemble those of compounds 1 and 3, suggesting that our  Coumarin is a well-known drug scaffold which displays versatile biological activities especially substantial anti-HIV activity (20). Notably, IN inhibitory activity has been reported for coumarin dimmers (1, Figure 2) (21) and 4-hydroxy-5-azacoumarin-3-carbox(thio)amides (2, Figure 2) (22). Furthermore, the anti-HIV activity of a bis-azo compound (FP-21399, Figure 2) has been evidenced (23). With regard to these data we selected coumarin scaffold as a core moiety to design new anti-HIV compounds. The decision for selecting novel designed compounds possess ability to be promising anti-HIV agents.

Materials
All reagents purchased from the Aldrich (USA) or Merck (Germany) Chemical Company and were used without further purifications.

General
Melting points (mp) were determined using a Thomas Hoover capillary apparatus (Philadelphia, USA). Infrared spectra were acquired on a Perkin-Elmer 1420 ratio recording spectrometer. A Bruker FT-500 MHz instrument (Brucker Biosciences, USA) was used to acquire 1 HNMR; Chloroform-D used as solvent. Coupling constant (J) values are estimated in hertz (Hz) and spin multiples are given as brs (broad singlet), s (singlet), d (double), t (triplet), q (quartet), m (multiplet), and br (broad). The mass spectral measurements were performed on an 6410 Agilent LCMS triple quadrupole mass spectrometer (LCMS) with an electrospray ionization (ESI) interface.

General procedure for preparation of 4-(5-(substituted thio)-1,3,4-oxadiazol-2-yl) aniline (8a-h)
To a solution of compound 7 (1.55 mmol) in methanol (4 mL) sodium hydroxide (10%, 0.12 mL) was added. The mixture was cooled in sodium nitrite aqueous solution at 0-5 °C. The final derivatives 11a-h were obtained through coupling reaction of the diazonium salts with 4-hydroxycoumarin (10) in aqueous sodium carbonate solution at neutral pH. an ice bath and substituted benzyl/alkyl halide (1.55 mmol) was added dropwise and stirred at room temperature. After completion of the reaction, the precipitate was filtered, rinsed with water, and crystallized from ethanol.
General procedure for the preparation of To a cooled solution of substituted compounds 8a-h (1 mmol) in 6 M hydrochloric acid (2.7 mL), 0.7 mL cooled solution of 10% sodium nitrite was added dropwise in order to keep the temperature less than 5 °C. Then, the mixture was poured into a solution of 4-hydroxy coumarin (10) (1 mmol) and sodium carbonate (0.05 g) in water (1.85 mL). The pH of the mixture was adjusted to 7 by addition of sodium acetate. The reaction mixture was stirred in an ice bath. After completion of the reaction, the precipitate was filtered, washed with water and crystallized from acetonitrile.

In vitro anti-HIV and cytotoxicity assays
Anti-HIV-1 activity of synthesized compounds was evaluated in single cycle replication assay which was evaluated in our laboratory and reported previously (32-35).

Biological evaluation
Compounds 11a-h were assayed for anti-HIV-1 activity in single cycle replication method which reflects only one round of infection. Anti-HIV activity of compounds was measured as inhibition rate (%) of p24 expression in Hela cells  11e and 11f). The shift of the chlorine atom from position 4 (compound 11f) to position 2 (compound 11e) was associated with decrease in anti-HIV activity (inhibition rate from 52 to 33%) and increase in cytotoxicity (cell viability from 100 to 75%). It can be concluded that 4-chlorobenzyl moiety at C-5 position of 1,3,4-oxadiazol ring was optimal in terms of anti-HIV activity and cell viability. cultures. Azidothymidine (AZT) was the positive control in the test. The cytotoxic properties were measured in the same cells by XTT proliferation assay and reported as percentage of cell viability. Results are summarized in Table 1. According to our results, all compounds indicated anti-HIV activity in the range of 5-79% at 100 μM concentration. AZT as positive control exhibited inhibition rate of 58% in our assay. Generally, all the tested compounds proved no considerable cytotoxic properties (cell viability > 74%) except compound 11a with 51% cell viability. Hence, it can be concluded that the anti-HIV-1 activity of the compounds was not a result of their cytotoxic effects. As it was expected, compounds containing halobenzyl substitutions exhibited better anti-HIV activity as compared to alkyl substitutions. Although compound 11a containing methyl substitution displayed inhibition rate of 79%, its cell viability was 51% which indicated that its anti-HIV-1 activity might be derived from cytotoxic property. Among alkyl substitutions, propyl (compound 11c) was reasonably well tolerated, giving similar results to the 2-chlorobenzyl substitution. Introduction of an unsubstituted benzyl group (compound 11d) was not well tolerated, resulting in reduction in anti-HIV-1 activity (inhibition rate = 6%). Introduction of fluorine atom in position 2 and 4 of benzyl ring (compound 11g and 11h, respectively) gave the similar inhibition rates (22 and 25%, respectively). However, compound 11g with fluorine atom in position 2 showed less %cell viability. The best anti-HIV-1 activity was observed with compounds containing chlorobenzyl substitution (compounds 11e and 11f). The shift of the chlorine atom from position 4 (compound 11f) to position 2 (compound 11e) was associated with decrease in anti-HIV activity (inhibition rate from 52 to 33%) and increase in cytotoxicity (cell viability from 100 to 75%). It can be concluded that 4-chlorobenzyl moiety at C-5 position of 1,3,4-oxadiazol ring was optimal in terms of anti-HIV activity and cell viability.

Molecular Modeling and Docking Study
A molecular docking study was conducted in order to evaluate synthesized compounds interactions with the integrase active site. To validate the docking study, the co-crystalized ligand, raltegravir was redocked under the same condition and superimposed on co-crystalized ligand pose (RMSD = 0.001), and displayed high affinity binding (-12.8 kcal/mol).
Docking results revealed that all the docked compounds accommodated well in the active site with affinity binding energy range from -9.0 to -9.8 kcal/mol (Table 1). Compound 11f showed high binding affinity (-9.8 kcal/mol). The 3D and 2D alignment of compound 11f in the active site was shown in Figure 3. The nitrogen atoms of 1,3,4-oxadiazole ring have been located in the 2.28 and 2.83 Å distance from Mg 2+ essential cofactors and was involved in metal acceptor interaction with the Mg 2+ ions. Phenyl and 1,3,4-oxadiazole motifs of the structure have been interacted via π-stacking bonding with deoxy adenosine (DA-17) presenting in the active site and could be expected to displace it. In addition, the 4-hydroxycoumarin motif was able to participate in π-stacking bonding with the Tyr212 residue.
Orientation of compound 11f and co-crystalized raltegravir in the active site shown in Figure 4 clearly emphasized that both compounds occupy the same place and are engaged in similar interactions with active site residues. Nitrogen atoms of 1,3,4-oxadiazole ring of compound 11f and chelating motif of raltegravir are nearly superimposable. 4-Chlorobenzyl group of compound 11f and 4-florobenzyl group of raltegravir have been similarly oriented towards the hydrophobic pocket. 4-Hydroxycoumarin motif of compound 11f and oxadiazole ring of raltegravir occupy the same area as well.
Collectively, the docking results of tested compounds were in a good agreement with IN inhibitors, suggesting that their anti-HIV activity may be due to IN inhibition.

In silico evaluation of drug-likeness
Oral bioavailability is considered as an important factor in design of new molecules as therapeutic agents. Thus, in drug design it is valuable to gain sufficient information about the molecular properties that limit oral bioavailability. There are some analyses that have correlated physical properties of molecules  with successful drug development. Among them, Lipinski′s rule of 5 (38) and Veber′s criteria (39) are the common principles utilized to evaluate the drug-like properties of a compound. In this study, Molinspiration online property calculation toolkit (40) was used to calculate the Lipinski′s molecular properties and the number of rotatable bonds (nROTB), together with the topographical polar surface area (TPSA; a sum of polar atoms′ surfaces: a descriptor for drug absorption, penetrability and bioavailability), the percentage of absorption (ABS%) calculated as (ABS% =109 -0.345 × TPSA) (41) and the molecular volume (a determinant of the transport characteristics). Lipinski parameters and Veber′s criteria of three FDA-approved INSTIs, raltegravir (RAL), elvitegravir (ELV), and dolutegravir (DTG) were calculated in a same way for comparison.
All compounds showed logP values ranged from 3.32-5.41, MW ranging from 380.38-490.93, HBA value of 8, HBD value of 1, nROTB values ranging from 3-7 (<10), TPSA values lower than 140 Å 2 (114.09 Å 2 ) which all suggesting that the evaluated compounds have acceptable flexibility and are expected to possess good permeability and oral bioavailability (39). According to these calculations, all the compounds displayed the percentage of absorption (ABS%) about 70 which is similar to ABS% of RAL, ELV, and DTG. Additionally, aqueous solubility (Log S) of the tested compounds was determined using OSIRIS Property Explorer (Data Warrior) software (42). As can be seen from Table 2, all the compounds exhibited weak to moderate aqueous solubility with Log S values ranging between -6.300 and -8.617 mol/L which are lower than solubility of RAL, ELV, and DTG.

Conclusions
The principal objective of this study was to investigate a new chemical scaffold with anti-HIV activity. A novel structure was constructed by hybridization of fragments from well-known anti-HIV-1 lead compounds. Key to the design was the introduction of a 1,3,4-oxadiazole ring as a chelating motif associated with a halobenzyl moiety to mimic pharmacophore of INSTIs. SAR analysis revealed that compound 11f containing 4-chlorobenzyl ring exhibited best anti-HIV-1 activity among synthesized compounds with inhibition rate of 52% near inhibition rate of AZT (58%). Docking studies showed that designed compounds bind to IN with a mode similar to co-crystalized raltegravir. Moreover, in silico computation of ADME properties of designed compounds showed non-violations of Lipinski′s rule of 5 and Veber′s criteria, suggesting that the designed compounds possess favorable physicochemical properties near to those of FDA-approved INSTIs, RAL, ELV, and DTG. However, the designed compounds displayed low aqueous solubility as compared to FDA-approved INSTIs which may be the main reason of their moderate anti-HIV-1 activity. So, these compounds could undergo further structural modifications to improve their anti-HIV-1 activity.